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ABSTRACT 

Recent observations of PSR J0045-7319, a radio pulsar in a close eccentric orbit with 
a massive main sequence B-star companion, indicate that the system's orbital period is 
decreasing on a timescale ~ 5x10^ years (Kaspi et al. 1996). Timing observations of PSR 
J0045-7319 also indicate that the B-star is rotating rapidly, perhaps close to its breakup 
rotation rate (Lai et al. 1995, Kaspi et al. 1996). For rapid (super-synchronous) prograde 
rotation of the B-star, tidal dissipation leads to an increasing orbital period for the binary 
system, while for retrograde rotation of any magnitude, the orbital period decreases with 
time. We show that if tidal effects are to account for the observed orbital decay of the 
PSR J0045-7319 binary, the B-star must have retrograde rotation. This implies that the 
supernova that produced the pulsar in this binary system likely had a dipole anisotropy. 

For a reasonably wide range of retrograde rotation rates, the energy in the dynamical 
tide of the B-star needs to be dissipated in about one orbital period in order to account for 
the observed orbital evolution time for the PSR J0045-7319 binary. We show, however, that 
the radiative dissipation of the dynamical tide in a rigidly rotating B-star is too inefficient 
by a factor of ~ 10'^, regardless of the magnitude of the rotation rate. We describe how, 
when the surface of the B-star is rotating nearly synchronously (which is expected from 
the work of Goldreich and Nicholson, 1989), the energy in the dynamical tide is dissipated 
in less than an orbital period, thus reconciling the theoretical and observed rates of orbital 
evolution. 

Nonlinear parametric decay of the equilibrium tide, for rigid retrograde rotation of 
the B-star, may also be able to explain the observed rate of orbital evolution, though the 
margin of instability is too small to draw definitive conclusions about the relevance of this 
process for the PSR J0045-7319 Binary. 
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§1. Introduction 

The recent discovery of two radio pulsars in binary systems with main sequence star 
companions has provided the opportunity to study tidal interactions with unprecedented 
precision (McConnell et al. 1991, Kaspi et al. 1994; Johnston et al. 1992). Of these two 
radio pulsars, the one in the Small Magellanic Cloud (SMC), PSR J0045-7319 (which has a 
spin period of 0.93 s), is particularly interesting for studying tidal interactions because the 
periastron separation between the neutron star and its companion, a main sequence B-star 
of nictss ^ 8.8Mq, is only about 4 times the B-star radius. The monitoring of pulse arrival 
times for the SMC pulsar has enabled accurate determination of several orbital parameters 
such as the period (51 days), eccentricity (0.808) and their time derivatives. Kaspi et al. 
(1996) report that the orbital period for the PSR J0045-7319 Binary is decreasing on a 
time scale of ~ 5 x 10^ years and that the orbital eccentricity is probably evolving more 
slowly. All relevant parameters for this system are included in table 1 for easy access. The 
work of Lai et al. (1995), Kaspi et al. (1996) and Bell et al. (1995) has also provided 
evidence that the B-star companion is rotating close to its break-up rotation rate and that 
its spin axis is misaligned with the angular momentum of the system. 

The evolution time for the SMC Binary's orbit due to the radiative dissipation of the 
dynamical tide (in a nonrotating B-star) is ~ 10^ years (see §2). It has been suggested (Lai, 
1996), however, that the observed rapid evolution of the SMC Binary can be accounted 
for by the radiative dissipation of the dynamical tide, provided that the B-star companion 
to PSR J0045-7319 has a rapid retrograde rotation rate of Cl^ ^ —0.4, where is 
the component of the stellar rotation rate perpendicular to the orbital plane in units of 
{GM^/Riy/"^. Retrograde rotation has the obvious advantage that the frequency of the 
tide seen from the rest frame of the B-star is considerably higher than for prograde rotation, 
causing excitation of very low order g-modes to large amplitudes. However, the dissipation 
times for low order g-modes, which set the rate at which orbital energy is extracted by the 
star, are considerably longer than the dissipation times of the higher order g-modes excited 
in a non-rotating star. This offsets the advantage of having increased mode amplitudes in 
a retrograde rotating star. For a concise and physically motivated derivation of this result 
please see Kumar & Quataert (1996). 

In the next section we describe the orbital evolution caused by the radiative damping 
of the dynamical tide in a uniformly rotating star, and discuss the probability that mode 
resonances may substantially speed up the rate of orbital evolution. In §3 we consider 
several alternate dissipation mechanisms, in particular nonlinear parametric instability 
(Kumar & Goodman 1996) and enhanced radiative dissipation in a differentially rotating 
star, in order to assess if they might be important for understanding the orbital evolution 
of the PSR J0045-7319 Binary. Please note that all frequencies in this paper, including 
rotational and orbital angular speeds, are given in micro-Hz (//Hz) unless explicitly stated 
otherwise. 

§2 Dynamical tide in a uniformly rotating star 
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Consider a binary system where the primary star has mass M*, and its companion, 
the secondary star, has mass M. The primary is assumed to be rotating as a sohd body 
with the component of its angular velocity normal to the orbital plane of fi* . 

The tidal forcing of a star predominantly excites gravity modes (g-modes) since they 
have periods comparable to the period of the tidal forcing. The g-modes in a uniformly 
rotating star are specified uniquely by three numbers: n the number of radial nodes, A, 
which is in general not an integer, is a generalization of the spherical harmonic degree £, 
and m, the azimuthal order. -"^ The mode amplitude in the rotating reference frame of the 
star, {a denotes the collective index n. A, m), for the dominant quadrupole tide, can 
be calculated by decomposing the tidal forcing function in terms of its Fourier coefficients 
and is given below (see Kumar et al. 1995, and Quataert et al. 1996 for details) 



^ V [^a ~ (k ~ |m|s)^]^ + d'^{k — |m|s)^ 



(1) 



where 
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Oo = a/ G{M + M^)/a^ is the orbital frequency, a is the semi-major axis of the orbit, 
D^^^^ are the fourier coefficients of the quadrupole tidal forcing function, is the mode's 
energy dissipation rate, and Qce is the overlap of the mode's normalized Eulerian density 
perturbation, Spa, with the quadrupole tidal potential. The negative sign in the exponent 
of equation (1) is considered when m is positive. 

We normalize our mode eigenfunctions so that the mode energy is equal to the square 
of the mode amplitude. Thus, the energy in mode a due to tidal forcing, as seen by an 
observer corotating with the star, is obtained by taking the square of the resonant term in 
the Fourier series, i.e., 



J a 



(3) 



[rl - {ka - \m\s)'^f + d1{ka - \m\s)'^ ' 

where ka. is a positive integer that minimizes the denominator, i.e., ka = [s\m\ ±ra\, 
where [x\ denotes the nearest integer to x. The net energy in the dynamical tide is 

It is easy to obtain D^^^^ numerically, and, for m = ±2, it peaks at A; ~ 2Q.p/Vto, 
where Vlp = (1 — e)~3/^(l-|-e)-'^/^Oo is the orbital angular speed of the star at periastron and 



^ Because the modification to the mode structure in the rotating star is not central to 
our analysis, we have relegated a discussion of it to Appendix B; in addition, we often use 
£ in place of A. 
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e is the orbital eccentricity. The peak value of -D^^^ is (1 + e) ^•\/l5/1287r(f2p/f2o) while 
for k > 5Qp/Qo, Df^^ ^ expi-l.Skflo/^p)/^^ - e)^/^ The m = Fourier coefficients, 
D^^\ are roughly constant for k up to 2Qp/Qo and decrease exponentially at larger k. 
Moreover, the D^^'^^ are typically smaller than the D^^^'^'^ by a factor of at least a few 
and so most of the tidal energy resides in the m = ±2 modes unless ~ Qp. The peak 

(22) 

of Df^ at /c ~ 20p/0o and the exponential fall off at higher harmonics implies that the 
mode with the largest energy has a frequency ~ ^ude (measured in a frame corotating 
with the star), where 0,tide = 2|rip — is twice the angular speed of the secondary at 
periastron, as seen in the rotating reference frame of the B-star. Higher frequency modes, 
while they have larger Qa, have less energy because Aa'' decreases exponentially for larger 
frequencies. Lower frequency modes have less energy because they have a larger number 
of radial nodes and thus smaller Qa- 

The effect of stellar rotation, for < f2* < ilp, is to decrease the energy in modes 
since the resonance is shifted to higher harmonics. The modes that are most efficiently 

(r) 

excited have uj^ ~ ^tide, and so E^-^^ decreases with increasing as lower frequency 
g-modes with smaller wavelengths and smaller Qa are excited. For fi* a; Qp or flude ~ 0, 
the energy in the dynamical tide reaches its minimum value. ^ 

For fl^ > flp, however, fltide increases with increasing O*, higher frequency modes 
are excited by the tidal forcing and so E^:^^^ increases. It is also clear from equation (3) 
that for retrograde rotation (il* < 0) the tidal frequency is larger and the energy in the 
dynamical tide is greater (as was pointed out by Lai, 1996). All of these effects are clearly 
seen in Figure (la), which is a plot of E^^j^, the sum of the energy in all modes (including 
m = & m = ±2), as seen by an observer corotating with the B-star, for the SMC binary 
system as a function of the B-star's rotation rate. 

The energy in a mode is proportional to 5r~^, where Sva = Va — \ka — \m\s\ (see eq. 
[3]) and so close resonances can significantly enhance the energy in the dynamical tide. 
The effect of resonances is evident in Figure (1). 

(r) 

The mode energy as seen by an observer corotating with the star, Ea , docs not 
include the work done by the tidal force on the velocity field associated with the rotation 
of the star, which changes the star's rotational energy. The mode energy as seen by an 
inertial observer includes this contribution and is given by Ea^ = E^^ + fi^Lo,, where 
Lq;, the angular momentum associated with a gravity wave of frequency cUq, and azimuthal 
order m can be shown to be equal to —mEa/<jJa- We note that E^ and Ua are frame 
dependent quantities, but the ratio Ea/uJa oc is frame independent. The net tidal 
energy as seen by an inertial observer, -E^^^g = 'Yla-^^^ ^ '^^ shown in Figure (lb) for the 
SMC binary system (we note that the results presented in fig. 1 were calculated using the 

^ This corresponds to the m = and m = ±2 modes being excited to comparable 
amplitudes. For Q.ude 7^ 0, the m = ±2 modes are excited to greater amplitudes than the 
m = modes. 
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entire fourier series expansion of E^^^^^ and E^'^^^, rather than just the resonant term; see 
Appendix A). 

Since (for m ^ 0) the azimuthal propagation of tidally excited waves, as seen by 
an inertial observer, is always in the direction of the orbital motion, the sign of m/u)a 

(i) 

is negative and is independent of fl^. This implies that Lq, and E^ have the same sign. 
Since the angular momentum deposited in the star is positive for ^ f^p, the tidal energy 
deposited in the star is also positive, which leads to a decrease in the orbital period with 
time. For > Op, on the other hand, the sign of the angular momentum deposited 
reverses and so there is a net transfer of energy from the star to the orbit. Physically, for 
ft* > ftp, the energy taken out of the spin of the star (which is being slowed down) exceeds 
the energy deposited in modes and so the orbital energy and period of the binary system 
increase with time. A more detailed discussion of these results is given in Appendix A. 

In order to calculate the orbital evolution we need to know the dissipation rate of the 
dynamical tide, which we calculate in the next section. 

2a. Radiative damping of g-modes 

The radiative damping times for the low order quadrupole g-modes of a 8.8 Mq, 
slightly evolved, main sequence star are given in Table 2 (the model was kindly provided 
by the Yale group). The dissipation was calculated by solving the fully non-adiabatic 
oscillation equations (see Unno et al. 1989 or Kumar 1994 for details). For n > 4, the 
mode dissipation time decreases rapidly with increasing mode order (decreasing frequency), 
roughly as a;~^-^. This scaling is explained below using a WKB analysis. 

The local radiative dissipation rate for the energy of a g-mode is (Unno et al. 1989) 

AT OAFr 

where A denotes a Lagrangian perturbation, T and Fr are the temperature and radiative 
flux, respectively, and e(r) is the mode's local energy density. Away from the turning 
points, this expression reduces to 

^ FRk^ f dlnT \ fdlnT \ 
^^^^ ~ pcldlnT/dr V dlnp )\ dlnp L ' ^ ' 



K^- [N' - col] W + 1)/-' - ^l/^'s] (6) 



where 

|V2r^/^ , 1^ /„2 , .2 /„2ll/2 

is the wave's radial wave number and Cg and are the soimd speed and Brunt- Vaisala fre- 
quency, respectively. The global radiative dissipation rate for a g mode (Fq,) is the integral 
of 7(r), weighted by e(r), between the lower and upper turning points: 

r^7(r)e(r), (7) 

J n 
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where r„, the outer turning point of the wave, is the radius at which uy^/cl{ru) = (-{i + 
l)/r^. Energy conservation implies that r^e{r)vg is independent of r (where Vg is the radial 
group velocity of the wave) and so the above equation reduces to 

Nu"^ r^pc^ [ J ' 

where N = J^" drN/r is the mean value of the Brunt- Vaisala frequency. Most of the 
contribution to the wave damping comes from a region near the upper turning point close 
to the surface of the star, which moves outward with decreasing frequency. Performing 
this integration for a polytrope of index 2 (a value appropriate for the outer envelope of 
the Yale group's B-star model we are using), we find that Fq, for g- modes scales as a; ~^ or 
(and also as £^). We note that for a fixed n, the value of Fq, for a quadrupole mode is 
independent of Jl*, cJq,, and m for a uniformly rotating star (which is why we have used £ 
instead of A in the above equations). 

2b. Standard dynamical tidal evolution for the PSR J0045-7319 Binary 

The rate of change of the orbital period for the SMC Binary can be calculated using 
the following equation 

Pnrh 'iEnrh 3 V — \„ 3 V— \ VrrEa^Ur'^ 



Porb '^Eqj-I) IE orb ^ "^Egrb ^ oj^^ — TTlVl 

where is the mode frequency in an inertial frame; the frequency in the rotating frame 
of the star is uj^'^ = uja^ — mQ^. The last equality in equation (9) follows from using 
the relation derived earlier between the mode energy in an inertial and rotating reference 

(r) 

frame. The mode energy in the rotating frame of the star, Ea , can be calculated using 
equation (3). The eigenfrequencies and the overlap of the modes with the tidal potential 
(Qa) are calculated using a 10 Mq B-star model of 5.3 solar radii (provided by the Yale 
group). The effect of the rotation of the B-star on the mode properties is included using 
the 'traditional approximation' (Bildsten et al. 1996; Chapman & Lindzen 1970; Unno 
et al. 1989; see Appendix B for details). The mass and the radius of the B-star are 
observationally estimated to be about 8.8Mq and 6i?0, respectively^ (Bell et al. 1995) and 



^ This radius is slightly smaller than that used by Lai (who took f» 6ARq), and so 
our tidal energies are smaller than his. It can be shown that the energy in the dynamical 
tide scales as Ej:2i^ oc Q^u;^ (see eq. [3]), where a is the mode that carries most of the 
tidal energy; the frequency of this mode in an inertial frame is ~ 2flp. Thus E^-^^^ oc 
(5a ~ {Rl/G){noi + 1)~^^, where (3 is equal to 1.5 (3) for a polytropic star of index 
3 (2), and 71^ is the number of radial nodes for the g-mode with the most energy. Since 
(n^ + l) oc {2np)-^y/GM^/Rl, we find that E^2le ^ RI'^^^/M^- thus for stars of polytropic 
index 3 (2) the tidal energy increases with the stellar radius (for fixed stellar mass) as i?^-^ 
{R\^); we have verified this dependence using numerical calculations. 



6 



so we multiply the mode frequencies and overlap integrals calculated using the Yale model 
by a factor of 0.78 and 1.36, respectively, to account for the difference between the stellar 
model and the observational parameters; note that Qq, oc i?^'^ and 00^ oc 

(M,/i?2)i/2. We 

have neglected the modification to the equilibrium structure of the star due to rotation in 
our calculations, which is likely to be a poor approximation since Vt^ for the SMC B-star 
is not much smaller than one. 

Figure (2) shows AEorb = PorbEorb and Porb/Porb for the SMC Binary as a function 
of the Bstar's rotation rate, fi* ; modes of m = ±2 as well as m = are included in these 
calculations.^ The results presented in this figure, as well as fig. 1, were obtained by 
artificially setting da = 0.2 in the calculation of the mode amplitude, thus limiting the 
strength of resonances; we will subsequently discuss the probability that closer resonances 
significantly enhance the energy in modes and thus decrease the orbital evolution timescale. 
The minima of E^f^^, -^tide ^Eorb (and thus the maxima of Porb/Porb) in Figures (1) 
and (2) correspond to the 'off-resonance' dynamical tide. They are equivalent to the 
calculation of the work done on the star in the absence of any oscillations, i.e., the result 
obtained by Press & Teukolsky 1975 (see also Lai 1996). 

For 17* < the tidal frequency is larger than for prograde rotation, and so the fre- 
quencies of the modes with the most energy are higher and their n value is lower. Thus 
the dissipation time of the dynamical tide is longer for retrograde rotation than for pro- 
grade rotation. For example, for fi* k, —0.3, the modes with the most energy are the 
gi-g^ modes, while for ~ 0, the modes with the most energy are the g^-gs modes. The 
dissipation times are longer in the former case by a factor of about 100 (see Table 2). This 
is why, although the energy in the dynamical tide is about 100 times larger for retrograde 
rotation than for no rotation (fig. 1), the energy dissipated, and thus the orbital evolution 
time, is comparable in the two cases. We note that 0,p/27r ~ 3.6//Hz, which accounts for 
the minimum in AEo^b and the maximum in Porb/Porb at Cl^ 0.2 (see fig. 2), since this 
is where flude ~ 0. We would also like to point out that the energy in the dynamical 
tide for is a few times 10"^*^ erg (fig. 1). If this energy is dissipated in an orbital 

period, the resulting orbital period evolution time is about 10^ years, which is close to the 
observed value. 

The observations of spin-orbit coupling in the timing data of PSR J0045-7319 provides 
evidence that the magnitude of the B-star's rotation rate is probably super-synchronous, 
that is, |ri*| > flp (Lai et al. 1995; Kaspi et al. 1996). Prograde rotation at il* > flp 
would lead to a transfer of energy from the spin of the B-star to the orbit, causing the 
orbital period to increase with time (see fig. 2 and the last part of the discussion following 
eq. [3]), which is inconsistent with the observations; prograde super-synchronous rotation 
of the B-star is thus ruled out. For il* ~ ^Ip (which is consistent with the apsidal motion 

^ We note that the results presented in Figures 1 & 2 were calculated using slightly 
more accurate expressions for Ea^ and Ea^ in which the entire fourier series, not just the 
resonant term, is kept (see Appendix A). 
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of the binary system) the energy in the dynamical tide is too small by about two orders 
of magnitude to explain the observed orbital evolution time, even if we assume that all of 
the tidal energy is dissipated in one orbit. These considerations force us to conclude that 
the rotation of the B-star must be retrograde. This is discussed further in §3. 

Our calculations indicate that the orbital evolution time due to the radiative dissipa- 
tion of the 'off-resonance' dynamical tide is at least three orders of magnitude too small 
to account for the observed Porb/Porb of the PSR J0045-7319 binary, regardless of the 
B-star's rotation rate. Is it possible, however, that the observed orbital decay rate of the 
PSR J0045-7319 Binary is because a low to moderate order g mode is highly resonant with 
the orbit {6ra ~ 10~^), thus increasing the energy in the dynamical tide by about three 
orders of magnitude and decreasing the orbital evolution timcscalc to about 5 x 10^ years? 
We think that this is highly unlikely. If the 'off-resonance' value of the energy in the dy- 
namical tide is Eor, then the probability that we are observing the system at a time when 
the energy in the dynamical tide is E is P{E) pa N^Eq^/E)^/'^, where N is the number of 
low to moderate order g modes which can be resonant with the orbit. This follows from 
comparing the time it takes the resonant mode to move off resonance with the time it takes 
the orbit to evolve between resonances. From Figures (1) and (2), E/Eor ~ 10^ in order 
to explain the observations and thus the probability that the orbital decay of the PSR 
J0045-7319 Binary is due to a resonant g mode is ~ 0.01%. This probability is nearly the 
same for = —0.3, 0. We thus conclude that it is highly unlikely that the orbital decay 
of the PSR J0045-7319 Binary is due to the radiative dissipation of the dynamical tide 
in a uniformly rotating star. We now consider several alternate dissipation mechanisms, 
including the effect of differential rotation on the damping of gravity waves. 

§3. Alternative damping mechanisms for tidal waves 

In the previous section we have argued that the linear radiative dissipation of the 
dynamical tide in a uniformly rotating B-star is incapable of explaining the observed orbital 
decay of the SMC Binary. The calculations of the energy in the dynamical tide are, in our 
opinion, relatively secure, but the understanding of dissipation mechanisms, which clearly 
affects our estimate of the orbital evolution time, requires closer scrutiny. Accordingly, 
we describe two dissipation mechanisms in this section. One of them, described below, is 
found to be extremely efficient in dissipating the energy of the dynamical tide when the 
star has some differential rotation. Another mechanism, nonlinear parametric coupling of 
the equilibrium tide to low frequency g-modes, is more likely to operate when the star is 
rigidly rotating. This is discussed in §3b. 

3a. The effect of differential rotation on the radiative damping of gravity waves 

The unusual (non-periodic) time residuals for the pulse arrival times of PSR J0045- 
7319 provides strong evidence for spin-orbit coupling in the PSR J0045-7319 Binary and 
suggests that the interior of the B-star is rapidly rotating with its spin axis inclined with 
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respect to the normal to the orbital plane (Kaspi et al. 1996; Lai et al. 1995). We 
can readily estimate the relative timescales for spin pseudo-synchronization and orbital 
circularization for the SMC system to see if the B-star should, from a theoretical point of 
view, be rotating pseudo-synchronously (i.e., with ^ ^p)-^ 

The rate of change of the stellar rotation due to the transfer of orbital angular mo- 
mentum to the star is 10*1/^^* ~ (Lorb/L^) {\Lorb\/ Lorb)i where Lorb/L^ ^ 40"-^. Using 
Lorb — Eorb/^p for the dynamical tide (see §2 and appendix A), as well as for the equi- 
librium tide^, we find that ^ 0.2 d~^{\Eorb\/Eorb)- Thus so long as the B-star 
is rotating near break up > 0.2), both the equilibrium and the dynamical tides lead 
to timescales for spin pseudo-synchronization which are comparable to the orbital circu- 
larization time. We therefore do not expect that the rotation rate in the interior of the 
B-star has changed appreciably, which is consistent with the observations. 

We do expect, however, that the dynamical tide should have forced the surface layers 
of the B-star into pseudo-synchronous rotation. The physical reason for this follows from 
the seminal work of Goldreich and Nicholson (1989), who showed that tidal waves deposit 
their angular momentum at the place in the star where they are dissipated. Since the 
dissipation rate is largest near the surface of the star (see §2b) and the moment of inertia 
of the surface region is a small fraction of the star's total moment of inertia, the surface 
of the star tends to be pseudo-synchronized very rapidly. For the SMC binary system 
with a rapidly retrograde rotating B-star (fi* pa —0.3), the modes with the most energy 
have frequencies ~ 17//Hz in the rotating reference frame. Figure (3) shows 7(r), the local 
radiative dissipation rate, for such a mode. The dissipation is concentrated in a layer 
~ 0.05-R* thick near the outer turning point, which occurs at r ~ 0.85-R*. The moment of 
inertia of this layer is ^ 10^ times smaller than that of the entire star and so we expect the 
rotation of the B-star at r 0.85i?* to be pseudo-synchronous. The dominant tidal waves 
are, however, evanescent at r > 0.85i?*, and lower frequency waves with outer turning 
points near the surface of the star (u;Q,/27r < S^uHz) carry insufficient angular momentum 
to bring the surface (r > 0.85i?*) into pseudo-synchronous rotation. The large positive 
entropy gradient in the radiative exterior of the B-star (the Brunt- Vaisala frequency is 
~ 70//Hz) has a strong stabilizing effect on shear instabilities which could otherwise retard 
the surface rotation and reduce the differential rotation of the region at r ^ 0.85i?*. 
Therefore, provided that magnetic stresses do not efficiently transport angular momentum 
between the stellar surface and the pseudo-synchronous layer, the rotation rate of the 
B-star at r > 0.85i?* is expected to be approximately the same as it was at the initial 
time. The optical linewidth measurement of Bell et al. (1995) gives a surface rotation 

^ The relative timescale calculation is significantly more secure than a calculation of 
the absolute timescale for pseudo-synchronization since it depends only on the ratio of the 
energy to the angular momentum in the tide. 

^ The equilibrium tide in the weak friction limit is discussed in detail in the seminal 
paper of Hut (1981), who also finds Eorb ~ ^pLorb- 
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frequency of ~ 4/ sin(z„s) /iRz, where ins is the angle between the spin axis and the hne 
of sight. If correct, this suggests that the observed surface of the star is perhaps not 
pseudo-synchronized. In what follows, we shall assume that there is a region below the 
surface of the star that rotates close to the pseudo-synchronous value and describe how 
this dramatically enhances the dissipation of the dynamical tide. 

We have calculated g-mode eigenfunctions in rigidly and differentially rotating stars 
and an example of the adiabatic energy flux for a g-mode is shown in Figure (4). The 
differential rotation was chosen to conform to that expected from the radiative dissipation 
of the dynamical tide in the B-star of the PSR J0045-7319 Binary, i.e., the thickness of 
the differentially rotating layer was taken to be O.li?* centered at r ^ 0.8i?*. We note 
that the energy flux is nearly constant across the differentially rotating layer, indicating 
that there is very little reflection of the wave at this layer; furthermore, the wavelength of 
the wave has decreased dramatically in the differentially rotating layer, which results in a 
strongly enhanced dissipation rate. 

The frequency of tidally excited gravity waves, as seen by an inertial observer, is ~ 2Qp 
(see §2a). If the component of the rotation rate of the star normal to the orbital plane a 
distance r from the center is ri*(r), then the frequency of the gravity wave in the local rest 
frame of the fluid is cUr ~ 2|f2*(r) — and its wavenumber is Q^^'^N{r) / {rujr). Thus, as 
the wave approaches the pseudo-synchronously rotating surface of the star its wavelength 
goes to zero and the energy it carries is entirely dissipated.^ Even if the rotation of the B- 
Star in the SMC Binary at r ^ O.SSi?* is not fully pseudo-synchronized, the dynamical tide 
will be completely absorbed as it approaches the surface provided that its local radiative 
dissipation rate becomes comparable to the wave frequency. We flnd that this occurs so 
long as the frequency of the gravity wave in the local rest frame of the star near r ^ 0.85i?* 
is less than about 2fj,liz, which requires that the stellar rotation frequency at r ^ 0.85i?* 
be within about l/iRz (or 30%) of its pseudo-synchronous value. This requirement, which 
is not particularly stringent, is weakened further if waves can travel closer to the surface of 
the star, such as when there is redistribution of angular momentum due to internal stresses 
in the outer envelope of the star. 

Let us assume that the energy input into the waves, inspite of the enhanced dissipation, 
is approximately equal to the energy in tidal oscillations calculated in §2a in the absence 
of resonances (the error due to this assumption is discussed below). From Figure 1 we 
see that this energy input is 2 x 10^^ ergs per orbit when the star is rotating with 
A* = —0.3 (fi*/27r — 6/iHz). Using Porb/Porb = —'^Eorb/^Eorb with an energy loss 

^ As the wave frequency goes to zero, the effect of the Coriolis force on the wavefunction 
becomes increasingly important and it might appear that this will adversely affect the 
dissipation of the wave. Using the 'traditional approximation' (Appendix B), however, it 
can be shown that, for a fixed wave frequency and a surface rotation rate which is slightly 
sub-synchronous, including the Coriolis force tends to increase the wavelength of the wave 
by only fa 20%, and thus the dissipation of the wave remains essentially the same. 
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per orbit of AEorb ~ Ell'^e ~ 2 X lO^^ ergs, we find an orbital period evolution time for 
the SMC binary of ^ 10^ years, which is a factor of about 5 smaller than the observed 
timescale. 

Thusfar we have ignored the modification to the g-mode eigenfunction in the differ- 
entially rotating star. As Figure (4) shows, in the outer part of the star where the stellar 
rotation rate approaches its pseudo-synchronous value, the wavelength of the gravity wave 
becomes very small. Thus, the coupling of gravity waves to the tidal forcing function is 
reduced, i.e., there is very little contribution to the overlap integral from the outer part 
of the star (r > 0.8i?*). The exact amount by which the overlap integral is reduced in 
the differentially rotating star depends on the details of the differential rotation and the 
modes which are excited, and is somewhat uncertain.^ Our estimates suggest that for 
the very low order modes excited in the rapidly retrograde rotating B-star, Qq decreases 
by a factor ~ 2 due to the differential rotation and so the energy in the dynamical tide 
decreases by about a factor ~ 4. This increases our estimate of the orbital period evolu- 
tion timescale of the PSR J0045-7319 Binary to ^ 4x10^ years, which is consistent with 
the observations. In fact the theoretically calculated orbital evolution times are consistent 
with the observations so long as 0.3 < J7* < 0. 

We note that prograde rotation of the B-star at il* ^ 0.4 gives an orbital evolution 
time of comparable magnitude, though with the opposite sign. Thus our reason for sug- 
gesting retrograde rotation is not to explain the short timescale for orbital evolution (which 
prograde rotation can also do), but rather to insure that the tidal interactions result in a 
decrease in the orbital period with time, as is found by the observations. For ~ 0.1 
(which is roughly the observational lower limit obtained from the apsidal motion of the 
system, cf. Kaspi et al. 1996) the energy in the dynamical tide is less than a few times 
Iq39 gj.g. ^j^g resulting orbital evolution time, assuming that all of this energy is dissipated 
in one orbit, is greater than about 3 x 10^ years, which is a factor of 50 larger than the 
observed value. For f» 0, on the other hand, the energy in the tide is about 2x10"^° 
erg and thus the orbital period is expected to decrease on a time scale of ~ 10^ years. 
This small rotation rate is, however, inconsistent with the observed apsidal motion of the 
system (Kaspi et al. 1996), which is why we are forced to consider retrograde rotation. 

It is straightforward to show that the rate of change of the orbital eccentricity is given 
by e ~ [(1 — e'^)/3e\{Porb/Porb)- Thus we expect ~ 3x10^ years, which is a factor of 
about 20 larger than the current lower limit of Kaspi et al. (1996). 



^ When the tidal frequency is much less than quadrupole f-mode frequency, almost all 
of the forcing by the tidal potential occurs near the interface of the convective core and 
the radiative exterior where the wavelength of the wave is the longest. Thus, differential 
rotation in the outer part of the star would not substantially modify the coupling of the 
wave to the tidal force. This limit is, however, not applicable to the very low order modes 
excited in the rapidly retrograde rotating B-star. 
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If magnetic stresses and/or instabilities are very efficient in redistributing angular 
momentum in the B-star so that, inspite of the tidally excited waves depositing angular 
momentum at r ~ 0.85i?*, the star continues to rotate almost rigidly, the dissipation 
of gravity waves will not be enhanced as described above. However, tides may still be 
efficiently dissipated by the nonlinear process discussed by Kumar and Goodman (1996). 
This is considered in some detail below. 

3b. Parametric coupling of the dynamical and equilibrium tides with g-modes 

When a primary mode, denoted by a and taken to be either the dynamical or the 
equilibrium tide, drives a mode of half its frequency (referred to as the daughter mode and 
denoted by /?) by nonlinear mode coupling, the process is called parametric instability. 
This process was proposed by Kumar & Goodman (1996) for damping the dynamical tide 
in late type stars, and was found to be much more efficient than conventional turbulent 
and/or radiative dissipation for close binary systems. We consider this process for the 
B-star of the SMC binary. 

The growth rate of the energy of the daughter mode, subject to parametric driving, 
is given by (Kumar & Goodman, 1996) 



J] 



(10) 



where E is the energy in the primary mode, Ao; = cJq: — 2a;/3, T/S is the linear energy dissipa- 
tion rate of the daughter mode, and Kq,/3^ is the nonlinear 3-mode coupling coefficient (the 
calculation of which is somewhat involved and is described in Kumar & Goodman). Thus, 
parametric instability sets in only if there are modes in the star which can simultaneously 
satisfy: 

|Aa;| < 'i\f2E^I'^Ko,s3fiU)fi (11) 

and 

< 6V^£;i/2««^^a;^. (12) 

These two conditions compete in the sense that higher i g-modes have smaller frequency 
spacings and are thus more likely to satisfy equation (11), but they have larger dissipation 
rates and are thus less likely to satisfy equation (12). 

The damping rate, F^, for g-modes of early type stars increases as f (see §2a), 
and so only low I daughter modes can satisfy equation (12). For retrograde rotation 
of = —0.3, the energy in the primary modes ((72 — (74) is i? ~ 2 x 10'^^ ergs (see Fig. 1), 
and the 3-mode coupling coefficient of these modes with daughters of half their frequency 
is «;a/3/3 ~ 10~^^ erg~^/^. Thus, we must find a daughter mode with |Aa;|/a;^ ~ 10~^ 
and < 3xl0~^ s~^ in order for the dynamical tide to be parametrically unstable. The 
small damping time required for the daughter modes implies that it must be of degree less 
than ~ 3. However, since the g-mode frequencies in the B-star of the SMC binary are 
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~ {n+l)~^^/I(T+Tj 27 //Hz, it is very unlikely (the probability is less than ^ 0.1%) that 
there is such a low degree daughter mode with a frequency within the required tolerance 
that can couple to the primary mode. Wc thus conclude that the dynamical tide in the 
PSR J0045-7319 Binary is stable against parametric instability. We show below, however, 
that the equilibrium tide is subject to the parametric instability. 

The equilibrium tide represents the hydrostatic response of the primary to the perturb- 
ing gravitational force of the secondary. The structure of the equilibrium tidal perturbation 
is similar to that of the / mode. The frequency of the equilibrium tide, ^ude — 21^^ — 0^,1, 
is, however, much smaller than that of the / mode. Thus the nonlinear coupling of the 
equilibrium tide to its daughter modes is analogous to the coupling of the / mode to very 
low frequency {up ~ ^ude <-^f) daughter modes, rather than to daughter's of half its 
frequency, as was the case for the dynamical tide. For Ufs ^ cu/, it can be shown that the 
coupling coefficient is proportional to the square of the daughter's normalized transverse 
displacement eigenfunction, which is proportional to (^^^. Thus, the coupling coefficient 
for the equilibrium tide is significantly larger than for the dynamical tide. For £ ^ 1 
daughter modes, ~ 4 x lO"^^^'^^ erg"-^/^, where up — Liip/27r is the daughter mode 

frequency in micro-Hz.^ 

The energy in the equilibrium tide near periastron is E^q k{R^/ Rp^riYGM'^ / ~ 
10^^ ergs, where RpeH ~ ^R* is the separation between the two stars at periastron and 
k ~ 0.008 is the apsidal motion constant of the B-star. The energy in the equilibrium tide 
falls off rapidly from periastron, so parametric instability can only set in and daughter 
modes can only grow near periastron. We thus define an orbit averaged growth rate for a 
daughter mode as = rjQo/^p: where rj is calculated using the energy in the equilibrium 
tide at periastron. If fj > Tp, the daughter mode has net growth over an orbit. 

For rapid prograde rotation of the B-star ($7* > Qp) the dissipation of the equilibrium 
tide causes the orbital period to increase with time, as was the case with the dynamical 
tide discussed in §2 (see appendix A and Hut 1981 for a detailed discussion). Thus, only 
nonlinear dissipation of the equilibrium tide in a rapidly retrograde rotating B-star can 
potentially explain the observed orbital period decrease of the SMC binary. 

For the SMC binary system, using the values given above for E^q and Kapp (and 
taking up ^ 8/uHz, appropriate for retrograde rotation at ~ —0.3) we estimate that 
7] ~ 10~^ s~^ and f) ~ 10~^ s~^. Thus we need to find daughter modes within about 0.2 
//Hz of 8 /tHz which is possible for £p > 2 when we consider the lifting of m-degeneracy 
by stellar rotation. We also find that ^ > so long as ip < 4 (using oc ij^ and the 
damping of quadrupole modes given in Table 2). Thus, we conclude that the equilibrium 
tide in the SMC Binary is subject to parametric instability which, as we see below, results 



^ For Up = 30/tHz (half the / mode frequency), this gives Kapp ~ 10~^^erg~^/^, which 
agrees well with the nonlinear coupling of the low order dynamical tide given above, as is 
expected since the equilibrium tidal perturbation is similar to that of the / mode. 
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in a rather efficient dissipation of its energy. 

The mean growth time of these daughter modes is fj ^ 100 days. The resulting 
nonhnear dissipation time for the equihbrium tide (tni) depends on the energy in the 
daughter modes, which is expected to be comparable to the energy in the equilibrium 
tide since the daughter modes are not likely to lose energy to granddaughter modes via 
parametric instability. Thus, tni ~ r/^^. The orbital evolution timescale which results from 
the nonlinear dissipation of the equilibrium tide is Porb/Porb ~ {Eorb/Eeq) (^p/^o) tni ^ 
77 (Eorb/Eeq) ~ 10^ ycars, in reasonable agreement with the observations. 

We feel, however, that the case for nonlinear damping in the SMC Binary, while 
encouraging, is not entirely robust. This is because the nonlinear growth rate of the 
daughter modes is comparable to their linear radiative damping rates and because there 
are not many daughter modes available for coupling with the equilibrium tide. Thus a 
change by a factor of a few in either or F/j, which is within the uncertainty of the 

stellar model, could modify our conclusion. 

4. Discussion 

The PSR J0045-7319 Binary in the SMC consists of a radio pulsar in a 51 day, highly 
eccentric (e=0.808), orbit with a B-star companion of mass ^ 8.8Mq. The orbital period 
of this system is observed to be decreasing on timescale of ~ 5x10^ years (Kaspi et al. 
1996), which is several orders of magnitude faster than that expected from the standard 
theory of the dynamical tide. Timing observations of PSR J0045-7319, as well as optical 
linewidth measurements, show that the B-star is rapidly rotating with a rotation speed 
that is perhaps greater than the orbital angular speed of the star at periastron (Lai et al. 
1995, Kaspi et al. 1996, Bell et al. 1995). 

We find that for rapid prograde rotation of the B-star (rotation frequency greater than 
the orbital angular speed of the star at periastron) the dissipation of both the equilibrium 
and dynamical tides lead to an increase in the orbital period of the binary system. This 
is because the tidal torque is slowing down the spin of the star and the resultant loss of 
rotational kinetic energy exceeds the energy deposited in the tide. Thus, we infer that the 
B-star must have retrograde rotation if tidal effects are to account for the observed period 
decrease of the PSR J0045-7319 Binary. Lai (1996) has suggested that rapid retrograde 
rotation of the B-star can explain the orbital evolution. We certainly agree with this 
conclusion. However, our reason for requiring retrograde rotation is entirely different from 
his. This point is clarified below. 

We find that for fi* = n^{Rl/GM^)^/'^ = -0.3 (0.0) the energy in the dynamical tide, 
as seen from an inertial frame, is about 2x10"^^ (2x10"^°) erg and the number of radial nodes 
for the most excited g-mode is 3 (6).^*^ The much larger energy in the dynamical tide for 

Our energies are smaller than those of Lai in part because we use different reference 
frames and stellar models, and also because we use a somewhat smaller stellar radius (see 
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rapid retrograde rotation was correctly pointed out by Lai. However, he took the damping 
time of g-modes to be independent of frequency (with a value of about 14 years) and 
concluded that the increase in the tidal energy for rapid retrograde rotation could account 
for the observed rapid evolution of the PSR J0045-7319 Binary. We have shown, however, 
that for a rigidly rotating star the mode dissipation time increases rapidly with frequency 
(roughly as w^) and so the high frequency modes excited in the retrograde rotating B-star 
are much less efficiently damped than the lower frequency modes excited in a non-rotating 
B-star. The damping time of the dynamical tide, in a rigidly rotating star, for = —0.3 
(0.0) is given in Table 2 to be about 1.1x10^(50) years (the dissipation was calculated by 
solving the fully nonadiabatic oscillation equations and the results are consistent with, for 
example, those of Saio & Cox 1980). Thus the orbital evolution time for rapid retrograde 
rotation of the B-star, inspite of the enhanced energy in the tide, is the same as for no 
rotation and is ~ 10^ years, or a factor of 10"^ longer than the observed value. We thus 
emphasize that retrograde rotation of the B-star is required, not because the standard 
theory for the dissipation of dynamical tide in a retrograde rotating star can explain the 
observed period decrease of the PSR J0045-7319 Binary (as Lai suggested), but rather 
because rapid prograde rotation of the B-star would lead to orbital period growth, which 
is inconsistent with the observations. In fact, with the enhanced dissipation mechanisms 
discussed below, the orbital evolution timescale for Cl^ ^ +0.4 is also consistent with the 
observations, but of course in this case the orbital period increases with time. 

The energy in the dynamical tide for Cl^ ?a —0.3 is about 10^^ erg (fig lb). If this 
energy is dissipated in about an orbital period then the resulting theoretical time scale 
for the evolution of the orbital period is similar to the observed orbital decay time for 
the SMC Binary. We have considered two different dissipation mechanisms that satisfy 
this property. One of these requires significant differential rotation of the B-star, which 
is expected when a star is not pseudo-synchronized. This is because the gravity waves 
deposit their angular momentum in the outer layers of the B-star where the dissipation 
occurs and where the moment of inertia is much smaller than in the interior (Goldreich 
& Nicholson 1989). For = —0.3, the dissipation of the dynamical tide in the B-star is 
concentrated at r 0.85/?* and so, provided that magnetic torques and/or instabilities do 
not force this region into solid body rotation with rest of the star, this region should be 
rotating pseudo-synchronously (i.e., with a rotation frequency close to the orbital angular 
frequency of the star at periastron, flp/27r 3.6/iHz). 

The frequency of the dynamical tide in an inertial reference frame is ~ 2Q,p and so the 
frequency in the local rest frame of the star is ~ 2[Qp~Q^{r)], where 0*(r) is the component 
of the stellar rotation rate at radius r normal to the orbital plane. Thus, the frequency and 
wavelength of the dynamical tide in the local rest frame of the star decrease as the wave 
approaches the psuedo-synchronously rotating layer. The decreasing wavelength results in 
the entire gravity wave energy flux being dissipated by radiative diffusion near the surface 

the footnote in §2b). 
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of the star, provided that the rotation frequency of the star at r ^ 0.85i2* is within about 
30% of the pseudo-synchronous value. This enhanced dissipation readily accounts for the 
observed orbital period decrease of the PSR J0045-7319 Binary so long as —0.3 < < 0. 
The resulting timescale for the evolution of the orbital eccentricity (1/e) of the PSR J0045- 
7319 Binary is ~ 3x10^ years, which is a factor of about 20 larger than the lower limit of 
Kaspi et al. (1996). 

The observed orbital period evolution can also be understood provided that the in- 
terior of the B-star is rotating very slowly (with the exterior close to the surface rotating 
pseudo-synchronously).^^ Slow rotation of the B-star's interior (|f2*| < 0.1) is, however, 
ruled out by the apsidal motion observations of Kaspi et al. (1996). Thus retrograde 
rotation of the B-star with Cl^ between f» -0.3 and -0.1 appears to be the only solution 
that fits all of the current observations of this binary system. 

In the likely scenario that the binary system was synchronized and circularized prior to 
the supernova that created the pulsar, the current retrograde rotation of the B-star implies 
that the supernova had a dipole asymmetry which imparted a net angular momentum kick 
to the system. This point was also made by Kaspi et al. on the basis of the misalignment 
of the spin and angular momentum axes. 

We note that in order for this mechanism to work the B-star must have a large 
gradient of differential rotation in the outer envelope, with the specific angular momentum 
decreasing outward in the star. This can lead to thermal instabilities such as the Goldreich- 
Schubert-Fricke instability, and it is unclear to us whether angular momentum will be 

mixed efficiently as a result (the star is, however, dynamically stable because of the large 
positive entropy gradient in the radiative exterior). This is an important issue which we 
have not addressed and which needs to be investigated by a careful nonlinear calculation. 

Another dissipation mechanism we have considered is nonlinear parametric instability 
of the equilibrium tide, which can operate when the differential rotation in the star is 
small. In this case energy is transferred, via resonant 3-mode coupling, from the large 
scale perturbation associated with the equilibrium tide to low degree g-modes of half the 
tidal frequency; the timescale for this process is found to be ~ 10^ days for the B-star of 
the SMC binary and the resulting orbital evolution timescale is a few hundred thousand 
years. The margin of instability in our calculation is, however, small and so we are not 
confident that this dissipation process is operating in the B-star of the SMC binary system. 
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which we found helpful to our thinking about this problem. We are very grateful to Jeremy 

Figure lb shows that the tidal energy for a slowly rotating B-star, Cl^ ^ 0, is about 
-|^q40 Qy-gg^ which is a factor of ~5 smaller than what is needed to explain the orbital 
evolution. The tidal energy, however, depends on the radius of the B-star, which is not 
precisely known; thus, if the radius of the star is some what larger than the value we have 
chosen, the tidal energy will be larger than the value shown in figure 1. 
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Appendix A 



In this appendix we calculate the relationship between the angular momentum and en- 
ergy deposited by the tidal force for the dynamical and equilibrium tides. The quadrupole 
tidal gravitational potential, as seen by an observer corotating with the primary star, is 
given by 



where 



and 



2m 



4. 

^ 5 °M + M* 



Y^*^{7T/2,(l)orb{t)-nj) 



iA2) 



(AS) 



The torque exerted on the primary star by the tidal force is only in the z direction (taking, 
for simplicity, the rotation axis to be perpendicular to the orbital plane) and is given by 



r X VU),5p, 



{AA) 



where 5p, the Eulerian perturbation to the density of the star, can be expanded in normal 
modes as 5p = AQ,5pQ,yg^(^, 0). Substituting the tidal gravitational potential (eq. 
[Al]) into equation (A4), we get 



Tz = -■i'P^'rnAaQaf2m- 



{A5) 



The non-dimensional function /2m can be expanded in a fourier series as (Kumar et al. 
1995; Quataert et al. 1996) 



/2m = exp{imQj) 



.k=l 



k=0 



{A6) 



l2±2\ 



where Im(D20) _ ^20 ^ ImiD^^^) =Re(C^±2) = 0, and, for A; > a few, Re(L'^ ; ~ 
— sign(m)Im(C^^^). Using equation (A6) one can expand the mode amplitude in a fourier 
series and obtain (c.f. Kumar et al. 1995) 



)2m ^ ■^2m^^Hknot - i0fe(m)] 

2m ,^2mNexp[i/cOot + i(/)fc(-m) 



r]k{m) 



+ {Dr-icr). 



Vki-m) 



{A7) 



where r]l{m) = [r^ — {k — ms)^]^ + d'^{k — ms)^ and the other symbols are the same 
as defined in equation (2) of the main text. Note that Aq, given above is as seen by an 
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observer corotating with the star. In what follows, we shall assume that ReiDl"^^) = 
— sign(m)Im(L)^^^) for all k. We find that this yields values for the energy and angular 
momentum deposited in the star which are accurate to within at least 5 %. Substituting 
the fourier series expansions of the mode amplitude and the tidal forcing function into the 
expression for the tidal torque [A5], we find that the angular momentum deposited in the 
star, in one orbit, is given by 

where we have summed over m = ±2 in arriving at the above expression. We note that 
m/oja has the same sign for m = ±2 modes (see eq. [A7] or eq. [2] of the main text) and 
so the modification to uja and Qa by the Coriolis force is the same for these modes (see 
Appendix B). 

The energy input rate into the star is given by 

^ = - J d^xpVU ■ V, {A9) 

where v, the fiuid velocity in the star, depends on the reference frame. In an inertial 
reference frame, Vi = Vr + x r, where Vr — ^^dAa/dt Yim{0,(j)) is the fluid 
velocity as seen by an observer corotating with the star (and is due to the tidally excited 
oscillations) and x r is the fluid velocity due to the rotation of the star (as seen by an 
inertial observer). It is straightforward to show from equations (A4) and (A9) that the 
energy input as seen by an inertial observer, E^l^^^, is related to the energy input as seen 

• ' (i) ' (r) 

by an observer corotating with the star [E^-^J by E^-^^^ = E^-^^ + Q^Tz. The change in 
the orbital energy per orbit includes both the energy input in the modes, -E^jje' 
changing rotational energy of the star (fi^r^) and is thus given by 

AEorb = - / dt 4^. (AW) 
Jo 

■ (r) 

We now calculate -E^jjg using the fourier series 

a a 

Substituting equations [A6] & [A 7] into the above equation we obtain 



^^tide = dt ii^ - Qn22^n22dn22 ^ * ~ ' 



+ 04 2^^n20'^n20"n20 2^r /^x.a'^ 



(A12) 
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where the first term in the above expression is due to the contribution of m = ±2 modes 
and the second (which is generally small compared to the first) is due to the m — modes. 
Note that the tidal energy as seen by an observer corotating with the star is positive 
definite (see Figure 1). The change in the orbital energy per orbit, AEorb = ~'^Etide ~ 
-(A£;([i + n,AL), is given by 

2^ Wn20^n20«^n20 2^ [^^^(^'^ 



For the results presented in this paper (e.g., Figures 1 &; 2), we have used the full 
fourier series expression for Efi^g, AE^rbj etc. (including maintaining the distinction be- 
tween and which, for simplicity, was dropped following equation [A7]). These 
expressions, however, contain summations over both mode order and fourier coefficients 
and are thus somewhat involved. We now discuss several approximations which allow for 
clearer physical understanding of the behavior of and relationship between AL, AE^^^^ 
and AEorb- 

The dominant contribution to the angular momentum and energy deposited comes 
from the resonant term in the fourier series, i.e., the integer k for which r]k{m) is mini- 
mized; we denote this value by k^- Keeping only this term in the fourier series, the above 
expressions for AL, AE^^^^ and AEorb simplify to 

AL ^ — ^ Ql22^l22dr.22^^^^{2Q^ - (^14) 



^^tide ^ Qn22^n22dn22J^I^-J-y^[^^U - ka) 

TT/?^ 2 2 T i^lD'^ 

+ 7^ 2^^n20^n200!n20|^^^-^fe, 

and 

AEorb = -A4i « ^ J]g^^2a;^22d„22i^^/c,(21], - /c,) 

^/5Vn2 2 , (Diy .2 
- 7)4- Z^^n20^n20dn20|^^^^-^/ea 



(^15) 



(^16) 



which are the expressions used in the main body of the text. We note that r]k{0) and ?7fc(2) 
will, of course, in general be minimized by different ka- These simplified expressions are 
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typically accurate to within 30 percent in calculating the energy and angular momentum 
deposited in the star. The only exception is when fip, at which point equations (A14)- 
(A16) can underestimate the energy and angular momentum deposited by up to a factor 
of a; 5. 

As a final simplification, we note that the modes with the most energy have Iml = 2 
and uja ~ 20p, where Ua is the mode frequency in an inertial reference frame. This is 
because, for k > Op/Og, the -D^^ are typically larger than the D'^^ by at least a factor of 
a few (> 3) and so the m = mode energies are in general negligible (with the possible 
exception of when ^ flp). Furthermore, D^^ peaks at A; 2Clp/Qo and decreases 
exponentially for larger k. Thus, modes with uja^ > 2Qp, while they have fewer radial 
nodes and hence larger Qq,, have less energy than modes with a;^*^ ^ 2Qp because of the 
rapid decrease of D^^ with k^- Similarly, for u>a^ < 2Qp, not only is D^^ smaller, but Qa 
is also smaller since the modes have more radial nodes. 

Thus the dominant contribution to AL and AEorb comes from the k^ ~ 2flp/Qo 
term in the mode summation, which dramatically simplifies the above expression. From 
equation (A14) wc sec that AL > for < Qp and AL < for > Op, indicating that 
the tidal force causes ^ ^ Qp (which is called pseudo-synchronous rotation). From 
equations [A14]-[A16] it follows that AE^ljie ~ OpAL and AE^l^^ (Op - n*)AL. These 
relations show that the energy in the dynamical tide is a frame dependent quantity. Finally, 
since AEorb = ~^Etide ~ — OpAL, we see that AEorb < for O* < Qp while AEorb > 
for O* > ftp (see Figure 2). The orbital period of a binary system thus increases due to 
the dissipation of the dynamical tide if > flp since the energy removed from the spin 
of the star (which is being slowed down by the tidal torque) exceeds the energy deposited 
in modes. 

We now show that the dissipation of the equilibrium tide leads to an analogous rela- 
tionship among E^-^^^, E^-^^, and r. For the equilibrium tide, we find it easiest to explicitly 
sum over the m = ±2 components in U, in which case 

U = -V{t) cos [2(0 - (l>orb{t) + nj)] + Vo{t), {A17) 

where V{t) ^ 2/3r'^Y^^2i^/2, 0)^2,2(6*, (p) / [R{t) / a]^ and the m = contribution to the tidal 
potential is given by Vo{t) = Pr'^Y^Q{7r/2,0)Y2,o{e, (l))/[R{t)/a\^. The equilibrium tide 
represents the hydrostatic response of the primary to the perturbing gravitational force of 
the secondary and so the time dependence of the equilibrium tide follows that of the tidal 
gravitational potential. Thus, the mode amplitude for the equilibrium tide, as seen by an 
observer corotating with the star, is given by (Kumar et al. 1995) 

Aa^ J (fxU{v,t)5p^. (A18) 

Substituting U and Aa. into equations (A4) and (A9) it is easy to show that, near peri- 
astron, e[^^^^ (Op — Q.^)tz and E^^^^^ OpT^, where we have taken 4'orb{t) ~ and 
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R{t) 0, which are vahd near periastron. For highly eccentric orbits, the energy in the 
equihbrium tide peaks strongly near periastron and so most of the contribution to AL and 
AEorb comes from near periastron. Thus, as with the dynamical tide, AEo^f) fa —QpAL 
for the equilibrium tide. 
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Appendix B 



Neglecting the rotational modification to the equilibrium structure of the star, the per- 
turbed mass and momentum equations, for a rotating star, in the Cowling approximation 
are 



and 



-u)'^p^ = -VSp + dpg + 2piu;^ x ft^, 



where 5 denotes an Eulerian perturbation and we have assumed a time dependence of 
exp(ia;Qt). The Centrifugal force does not appear in equation (B2) because it has no Eule- 
rian perturbation.^^ The components of the momentum equation for adiabatic oscillations 
can be written as 



-U = 



d6p g6p 



dr 



N'^^rP + 2ia;pn*^<^ sin 9, 



2 ^ 1 dSp „ . ^ ^ ^ 
-oj pt,e = -z~J7r + 2za;pS2*^<^cos6/, 



r d9 



and 



1 dSp 



2iujpVl^ {^r sin 9 -\- cos 6*) . 



(S3) 



(S4) 



(S5) 



rsin^ (i0 

where we have taken = f2*z. In the region of wave propagation ~ N^r/^i and thus 
the buoyancy term in the radial momentum equation is larger than the Coriolis term by a 
factor of N /Q,^.-. we thus neglect the latter. Away from the equator, sin0 ~ {^e cos 9)uj/N 
and thus we neglect the radial displacement in the component of the momentum equation. 
These approximations constitute the 'traditional approximation' (Bildsten et al. 1996; 
Chapman & Lindzen 1970; Unno et al. 1989) and are only valid when ojo, ^ N and 

Solving for and in terms of Sp and taking the perturbed variables to be propor- 
tional to exp(zm0) we obtain the following coupled radial equations 



dr ^ pCg 



1 



^ = 

ci 



(56) 



and 



dSp gSp 
dr 



p(a;2 - N^)^r = 0, 



{B7) 



■^^ Equivalently, the Centrifugal force evaluated at the perturbed position of the fluid is 
identically canceled by the equilibrium structure of the star evaluated at the perturbed 
position. 
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where the angular eigenfunctions (called the Hough functions) and eigenvalues (A) are the 
solutions to the following eigenvalue problem 



LqH — — XH, 



{B8) 



where 



— d I 1 — x'^ d 




(S9) 



{1-x 



^)(1 — g^x^) dx \_1 — q'^x'^ dx 



— qm 



(1 - q'^x^) 



2 ' 



where x = cos 9 and where the angular operator (and thus A and the Hough functions) 
depends on the stellar rotation rate and mode frequency through the parameter q = 
2n^/uJcf We note that, for a fixed |m|, \u!a\, and it is the sign of mQ^/ua that 

determines the effect of the Coriolis force. 

Figure (5a) shows A as a function of g = 2fl^/iVce for the even Hough function with 
the smallest value of A, which is a generalization of £ = 2. For m = —2, A rn? = 4 
for q^ 1 and thus the mode structure of the m = —2 modes is only weakly modified by 
the rotation of the star. More detailed discussion of the Hough functions is given in the 
above references. In what follows, we focus on the aspects of the modified mode structure 
relevant to tidal excitation. 

The overlap integral for a mode can be written as 



It can easily be shown using a WKB analysis, or verified numerically, that Qo,,r-, the radial 
contribution to the overlap integral, is only a function of the number of radial nodes of 
the mode, i.e., for a fixed n, Qa,r is independent of Q.^, uj^, A and m. The overlap of the 
even Hough function with the smallest value of A (the generalization of £ = 2) with the 
quadrupole spherical harmonic, Q angle-, is shown in Figure (5b) as a function of q. For 
g < 1, the associated Legendre polynomials are very good approximations to the Hough 
functions, but for g > 1 the Hough functions are concentrated near the equator. 

Using the information in Figure (5), there is a simple prescription for calculating g- 
mode frequencies and overlap integrals of a rotating star in the 'traditional approximation'. 
To calculate the frequency of a quadrupole g^ mode, take the value of q^ = 20* /o;^ 
(superscript o denotes the value for non-rotating star) and determine A from Figure (5a). 
The frequency in the rotating star is, to first approximation, given by cu^A^/^/v^. This 
frequency can be used to re-evaluate qa and in turn A and Wq, and this process can be 
iterated to the desired accuracy. The angular overlap integral, Q angle-, can be determined 
from Figure (5b) using the resulting value of qa, and, since Qa,r is a function of n alone, 
independent of O*, we obtain for the g-modes of the rotating star. 




(SIO) 
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Table 1 

Parameters for the SMC binary system 



Orbital period {Porb) 


4.421 X 10^ s (51.17 days) 


e 


0.8079 


Porb/ Porb 


-4.63 X 10^ yr. 


l/|e| 


> 1.5 X 10^ yr. 


Op/27r 


3.61 /uHz 


M, 


^ 8.8 Mq 




^ 6.0 Rq 




10.1 //Hz for = 0.5 
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Table 2 



Frequencies and radiative dissipation times (energy) for the 
low order quadrupole g-modes of a non-rotating B-star 
of mass 8.8 Mq and radius 6.0 Rq 



n (radial nodes) 


frequency (//Hz) 


(yrs) 


1 


34.66 


8.7x10^ 


2 


23.04 


2.5x10^ 


3 


17.17 


1.1x10^ 


4 


12.58 


4.6x10^ 


5 


11.21 


1.6x10^ 


6 


9.58 


47.9 


7 


8.43 


14.5 


8 


7.56 


5.6 


9 


6.83 


2.5 


10 


6.23 


1.1 
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Figure Captions 



FIG. 1. — (a) The energy in the dynamical tide for the B-star of the SMC Binary, as seen 
by an observer corotating with the B-star (E^ll^), as a function of the B-star's rotation rate 
(in units of ^GM^/Rl). (b) The energy in the dynamical tide for the B-star of the SMC 
Binary, as seen by an inertial observer {E^^i^^). Note that -E'fj'^g > for all rotation rates 
while E^^^^ > for < J^p, i.e., < 0.2 (soHd Hue) but E^^i^ < for Q,. > Qp (dashed 
line). All relevant £ — 2 g-modes of m = 0, ±2 were included in the calculations. These 
calculations (as well as those of Fig. 2) were made using the full fourier series expressions 
for E^ll^^ and E^^^^ given in Appendix A. 

FIG. 2. — (a) The change in the orbital energy per orbit, AEorb, and (b) the resulting 
orbital evolution timescale, Porb/Porb, for the radiative dissipation of the dynamical tide 
in a rigidly rotating B-star of the PSR J0045-7319 Binary as a function of the B-star's 
rotation rate (in units of ^^/GM^JE^); the calculation included all quadrupole g-modes of 
m = ±2 as well as m = 0. For < Op, i.e., < 0.2, AEort and Port/ Port are negative 
(solid line) while for > Op, AEorb and Porb/Porb are positive (dashed line); please see 
the text for a physical explanation of this result. 

FIG. 3. — The local radiative dissipation rate for a quadrupole g-mode of frequency ~ 
ITfiRz for the B-star of the SMC Binary; this is the mode predominantly excited for 
rapid retrograde rotation of the B-star (O* = —0.3). The dissipation peaks near the outer 
turning point of the mode, which occurs at r ~ 0.85il*, and so the angular momentum 
and energy of the mode are deposited well beneath the surface of the star. 

FIG. 4. — The adiabatic energy flux in a g-mode of the B-star in the SMC Binary for rigidly 
(dashed line) and differentially rotating (solid line) stars which have the same rotation 
frequency in the interior. The differential rotation is centered at r ^ O.Si?*, occurs over 
a width of ~ O.li?*, and is such that the minimum wave frequency (which is ~ 17/uHz in 
the core) is ~ 1/xHz at r ~ 0.8-R*. The small wavelength in the differentially rotating layer 
implies that the mode dissipation is strongly enhanced. Note that energy fliix is nearly 
constant across the differentially rotating layer, or in other words the differential rotation 
does not significantly reflect the wave. 

FIG. 5. — (a) The smallest value of the angular eigenvalue A as a function oi q = 2fi,^/(jja 
for even parity modes (these are the generalization of the i = 2 modes for a rotating star 
and so A — s> £{£ + 1) = 6 as g — * 0). (b) Q angles the overlap of the angular eigenfunction 
in the rotating star with the quadrupole spherical harmonic, for the same modes as part 
(a). Since the radial contribution to the overlap integral is, for a fixed number of radial 
nodes, independent of O*, Q angle give the variation of the overlap integral with the stellar 
rotation rate. 
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Figure 1 
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